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INTRODUCTION 


1  . 

Although  their  development  is  comparatively  recent,  super 
conductors  have  already  been  applied  to  power  generation, 
computers,  and  metrology.  Additional  commercially  promising 
applications,  as  well  as  their  possible  use  in  military  systems, 
are  currently  being  explored.  However,  to  make  superconductors 
more  attractive  to  potential  users,  cryogenic  refrigerators  need 
to  be  developed  that  are  capable  of  generating  (efficiently  and 
reliably)  the  required  ultra-low  temperatures. 

One  of  the  most  promising  superconductor  applications  is  the 
detection  of  magnetic  fields.  However,  it  is  this  application 
where  the  need  for  a  fieldable  cryogenic  refrigerator  is  perhaps 
the  most  acute,  and  where  the  superconductor/refrigerator 
interface  is  the  most  demanding.  Under  contract  to  the  Office 
of  Naval  Research  (ONR),  Philips  Laboratories  studied  the 
applicability  of  the  Stirling  refrigeration  process  to  the 
cooling  of  magnetometers,  and  has  generated  a  preliminary 
refrigerator  design  capable  of  meeting  ONR's  specifications. 

The  study  had  two  major  aspects:  thermodynamic  and  interactive. 
The  first  addressed  the  thermal  and  physical  characteristics  of 
the  refrigerator.  The  goal  was  to  determine  whether  the  cooling 
required  by  a  magnetometer,  specified  as  50  mW  at  10*K,  can  be 
generated  by  a  refrigerator  having  an  efficiency,  weight,  and 
size  compatible  with  military  applications,  i.e.,  have  relative¬ 
ly  low  power  input,  be  comparatively  light,  and  be  small. 

The  second  aspect  of  the  study  dealt  with  those  facets  of  re¬ 
frigerator  design  which  are  thought  to  affect  the  performance 
of  the  device  to  be  cooled,  viz.,  the  magnetometer.  The  goal 
was  to  assess  the  nature  of  the  magnetometer/refrigerator 
interface,  i.e.,  the  cold  surface  motion,  the  stability  of  the 
output  temperature,  and  the  magnitude  of  the  magnetic  fields 
generated  by  the  refrigerator  at  its  cold  region.  More  specif¬ 
ically,  the  objective  was  to  determine  whether  these  factors 
would  be  compatible  with  the  desired  magnetometer  output. 
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The  thermodynamic  aspects  of  the  study  had  no  significant 
uncertainties.  Existing  analytical  tools,  validated  in  labora¬ 
tory  prototypes  and  in  commercially  available  hardware,  were 
applied  to  a  combination  of  proven  concepts  and  novel  approaches 
to  yield  design  data  which  was  then  translated  into  a  layout 
drawing.  Specifically,  the  Philips  Stirling  Computer  Program 
was  used  to  optimize  the  critical  parameters  of  a  linearly 
driven,  triple-expansion  configuration.  Where  uncertainties  did 
exist,  the  assumptions  that  had  to  be  made  were  "conservative", 
to  provide  confidence  that  the  performance  predicted  by  the 
study  could  be  attained  in  working  hardware. 

The  interactive  aspects  of  the  study  were  less  precise  for  three 
major  reasons:  the  exact  geometry  and  the  physical  character¬ 
istics  of  the  magnetometer  to  be  cooled  were  not  inputs  to  the 
study;  prior  experience  with,  and  current  knowledge  in,  magnetic 
field  generation  in  cryogenic  hardware  are  not  sufficient  for 
rigorous  quantitative  work;  and  the  magnitude  of  all  three 
parameters  of  interest,  i.e.,  vibrations,  temperature  varia¬ 
tions,  and  induced  magnetic  fields,  are  subject  to  parametric 
trade-offs.  Such  a  parametric  study  can  best  be  done  as  part  of 
the  detailed  design,  i.e.,  in  the  next  phase  of  the  program. 
The  nature  of  that  phase  of  the  program  was  the  last  item 
addressed  by  the  study. 

This  report  outlines  the  results  of  Philips  Laboratories' 
analysis. 
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2. 


TECHNICAL  DISCUSSION 


2 . 1  Design  Concept 

A  Stirling  refrigerator  generates  cold  by  periodically  expanding 
and  compressing  an  ideal  gas.  The  expansion  takes  place  at  the 
desired  low  temperature.  The  compression  is  maintained  at 
ambient  temperature  with  the  aid  of  an  external  coolant.  The 
working  gas  is  usually  helium. 

The  expansion  and  compression  of  the  helium  is  accomplished  by 
the  reciprocating  motion  of  a  piston,  complemented  by  the 
out-of-phase  reciprocation  of  a  displacer.  The  latter,  as  its 
designation  implies,  displaces  gas  from  the  cold  region  to  the 
ambient  temperature  region  of  the  machine  (and  vice  versa) 
through  a  thermal  regenerator. 

The  efficiency  of  a  Stirling  refrigerator  depends  to  a  very 
large  extent  on  the  manner  in  which  the  required  expansion, 
compression,  and  regeneration  is  accomplished.  For  the  ONR 
application.  Philips  Laboratories  selected  a  triple-expan¬ 
sion  configuration  which  has  three  expansion  stages  and  three 
associated  regenerators.  The  three  expansions  are  produced  by  a 
single  displacer  (having  three  steps)  which  contains  the  three 
regenerators . 

The  efficiency  of  a  cryogenic  refrigerator  also  depends  on 
the  effectiveness  with  which  the  compression  heat  is  dissipated. 
In  the  proposed  design,  that  heat  will  be  transferred  to  a 
liquid  coolant  in  a  counterflow  heat  exchanger.  The  proposed 
heat-dissipation  method  has  been  used  often  in  the  past. 

The  physical  and  operational  chracteristic  of  a  Stirling  re¬ 
frigerator,  viz.,  it's  size,  weight,  endurance,  and  reliability, 
depend  to  a  large  extent  on  the  manner  in  which  the  piston  and 
displacer  are  driven.  The  proposed  design  uses  a  linear  drive 
in  which  there  is  no  need  to  translate  rotary  motion  into 
rectilinear  motion.  Because  its  motion  is  purely  linear,  the 


proposed  drive  requires  no  bearings.  Also,  the  seals  on  the 
piston  and  displacer  are  not  subjected  to  the  side  forces 
normally  present  in  a  rotary/  linear  drive  approach. 

A  critical  requirement  in  the  application  at  hand  is  that  the 
refrigerator/magnetometer  interface  be  thermally  efficient.  The 
most  difficult  aspect  of  this  interface  is  the  motion  (or 
vibration)  of  the  cold  surface,  induced  by  the  out-of-phase 
reciprocation  of  the  moving  masses.  To  eliminate  such  motion, 
or  reduce  it  to  a  negligible  level,  the  proposed  refrigerator 
design  will  be  passively  balanced  by  a  mass  equivalent  to,  and 
moving  in  opposition  to,  the  piston/displacer  combination. 

Many  features  of  the  proposed  design,  such  as  the  triple-expan¬ 
sion  configuration,  the  linear  drive,  and  the  balancing  ap¬ 
proach,  have  been  demonstrated  in  operational  hardware.  The 
characteristics  of  these  features  are  discussed  in  detail  in  the 
three  sections  which  follow.  Two  other  items,  however,  the 
coldest  regenerator  and  the  piston  seal,  although  operationally 
demonstrated  in  the  past,  are  thought  to  have  a  significant 
influence  on  the  reliability  of  the  proposed  design,  and  are 
therefore  discussed  in  a  section  devoted  to  the  critical 
aspects  of  that  design. 

2.1.1  Triple-Expansion  Configuration 

To  generate  the  required  50  mw  of  refrigeration  at  10’K  with  a 
relatively  high  efficiency,  the  proposed  design  is  a  triple¬ 
expansion  configuration.  A  detailed  description  of  multiple 
expansion  configurations  is  given  in  Appendix  A;  Philips 
Laboratories'  specific  experience  with  triple-expansion  designs 
is  also  presented  in  this  appendix.  In  principle,  the  triple¬ 
expansion  process  can  be  extrapolated  one  step  further  (i.e.,  to 
four  expansions)  to  attain  temperatures  below  8.5*K,  a  possible 
future  requirement. 


Figure  1  depicts  the  proposed  triple-expansion  approach  by 
showing  a  cross-section  of  the  cold  side,  i.e.,  the  combination 
of  cold-finger  (schematically),  displacer,  and  three  (3)  re¬ 
generators.  The  cold  finger  is  shown  in  greater  detail  in  the 
layout  drawing  of  Par.  2.1.5.  The  manner  in  which  the  displace’* 
is  integrated  into  the  rest  of  the  refrigerator  is  shown  on  that 
same  drawing. 
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Figure  1.  Proposed  triple¬ 
expansion  configuration. 


As  the  arrows  in  Figure  1  indicate,  the  helium  working  gas  flows 
from  the  compression  to  the  expansion  space  (and  back)  through 
three  regenerator  matrices.  At  three  points  on  this  path,  the 
gas  is  expanded.  It  is  during  these  expansion  processes  that 
refrigeration  is  produced . 


The  first  (warmest)  regenerator  matrix  is  made  of  phosphor 
bronze  mesh  with  a  filling  factor  of  30%.  It  is  31  mm  in 
diameter  by  49  mm  long.  Its  cold  end  temperature  is  100®K. 


The  second  (intermediate)  regenerator  matrix  is  also  made  of 
phosphor  bronze,  but  with  a  filling  factor  of  42%.  It  is  21  mm 
in  diameter  by  27  mm  long.  Its  cold  end  temperature  is  50*K. 

The  third  and  coldest  regenerator  matrix  is  made  of  lead 
spheres  which  are  140  in  diameter.  The  matrix  is  14  mm  in 
diameter  by  46  mm  long.  Packed  spheres  (of  any  size)  yield  a 
63%  (nominal)  filling  factor.  The  configuration  of  the  lead 
(i.e.,  whether  spheres  or  another  geometry)  in  the  matrix  is 
critical  to  the  life  and  reliability  of  the  refrigerator.  This 
aspect  of  low  temperature  regeneration  is  discussed  in  detail  in 
Par.  2.2.1. 

The  helium  working  gas  is  prevented  from  flowing  through  the  gap 
between  the  outside  of  the  displacer  and  the  cold  finger,  and 
thus  from  short-circuiting  the  regenerators,  by  the  three 
sleeve-type  seals  shown  in  Figure  1.  Seal  effectiveness  and 
endurance  are  discussed  in  Par.  2.2.2. 

The  analysis  which  led  to  the  design  shown  in  Figure  1  was 
based  on  the  requirement  for  50  mW  (minimum)  of  refrigeration  at 
10”K,  and  zero  refrigeration  at  8.5“K.  No  net  refrigeration  is 
available  at  the  intermediate  temperatures  (i.e.,  at  50°K  and 
IOO'K).  However,  these  intermediate  temperature  levels  could  be 
used  to  cool  radiation  shields  which  would  reduce  the  load  on 
the  coldest  (i.e.,  lO'K)  region.  However,  an  accurate  trade-off 
study  of  such  a  design  can  be  rigorously  addressed  only  when  the 
geometry  of  the  item(s)  to  be  cooled  is  precisely  known. 

It  should  be  noted  that  the  displacer/regenerator  subassembly 
reciprocates  500  times  per  minute  at  a  maximum  amplitude  of  7  mm 
within  the  cold  finger.  Also,  under  the  stresses  imposed  by  the 
varying  pressure  in  the  refrigerator,  the  cold-finger  wall 
strains  cyclically.  Both  the  intended  and  unintended  motions 
(i.e.,  the  displacer  movement  and  the  cold-finger  strain, 
respectively)  are  expected  to  generate  undesirable  magnetic 
fields,  as  discussed  in  Par.  2.3. 
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In  summary,  the  combination  of  the  cold- side  geometry  depicted 
in  Figure  1  and  the  matrix  materials  selected  for  the  three 
regenerators  is  expected  to  produce  the  required  refrigeration 
at  the  specified  temperatures.  However,  the  nature  of  the  lead 
matrix  and  the  choice  of  cold-side  structural  materials  require 
further  study,  as  discussed  in  Pars.  2.2.1  and  2.3. 

2.1.2  Linear  Drive 

To  impart  the  required  rectilinear  motion  to  the  piston  and 
displacer,  the  proposed  design  uses  a  purely  linear  drive, 
i.e.,  one  devoid  of  rotary-to-linear  conversion.  The  first  such 
drive  was  designed  at  the  Philips  Research  Laboratories  in 
Eindhoven,  The  Netherlands,  and  applied  to  a  Stirling  re¬ 
frigerator  (designated  Model  MC-80 )  which  produces  1  w  of  cold 
at  n°K.  Currently,  several  linearly  driven  refrigerators  with 
linear  drive  configurations  are  being  developed  within  Philips, 
including  two  at  Philips  Laboratories  in  Briarcliff  Manor. 

A  linearly  driven  refrigerator  is  often  referred  to  as  having  a 
free  displacer.  Indeed,  the  displacer  is  actuated  by  the 
pressure  difference  created  as  a  result  of  the  flow  losses 
through  the  regenerator  (rather  than  mechanically).  The  gas 
flow  is,  of  course,  induced  by  the  motion  of  the  piston  which  is 
driven  by  a  linear  electric  motor. 

A  schematic  drawing  of  the  proposed  linear  electromagnetic  drive 
system  is  given  in  Figure  2.  The  major  element  of  the  system, 
viz.,  the  motor,  consists  of  a  moving  coil  which  is  directly 
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Figure  2.  Linear  electro¬ 
magnetic  drive. 


coupled  to  the  piston  and  free  to  reciprocate  in  a  gap  formed  by 
a  stationary  permanent  magnet  and  an  iron  yoke  assembly.  In 
such  a  configuration  the  gap  length  to  gap  width  ratio  is  large, 
which  permits  a  nearly  uniform  magnetic  field  to  be  established 
across  the  gap.  Within  that  gap,  a  coil  can  reciprocate  as  the 
result  of  an  applied  ac  current.  The  piston  and  linear  motor 
are  both  sealed  in  a  hermetic  envelope,  with  the  electric  power 
transmitted  through  electric  feedthroughs.  The  center  position 
of  the  coil-piston  assembly  is  held  and  controlled  by  a  small 
mechanical  spring. 

It  should  be  noted  that  linear  motors  can  also  be  designed  with 
moving  magnets  (rather  than  moving  coils),  in  which  case  the 
need  for  flexing  electrical  leads  is  eliminated.  However, 
moving-coil  type  motors  are  more  efficient,  a  parameter  of 
importance  in  this  application.  The  motor  proposed  for  this 
application  has  an  efficiency  of  70%  at  maximum  amplitude.  The 
flexing  leads  are  not  expected  to  affect  the  life  and  reli¬ 
ability  of  the  proposed  design;  experience  with  the  MC-80,  as 
well  as  life  test  results  with  othc-i:  units,  suggests  a  flexing- 
lead  endurance  life  considerably  in  excess  of  those  expected  in 
this  application. 

As  indicated  above,  the  displacer  is  "free",  in  as  much  as  its 
motion  is  not  mechanically  imposed.  To  restrict  its  travel  and 
to  store  energy  as  a  means  of  enhancing  refrigerator  efficiency, 
the  ambient  temperature  extremity  of  the  displacer  is  anchored 
to  the  housing  through  a  spring,  as  shown  in  the  layout  drawing 
given  in  Par.  2.1.5. 

It  should  be  noted  that  there  is  a  significant  difference 
between  a  conventional  refrigerator  drive  and  the  linear  ap¬ 
proach  proposed.  In  the  former,  the  dynamics  of  the  system  are 
set  by  the  parameters  of  the  mechanical  drive;  therefore,  the 
optimization  of  such  a  system  deals  with  thermodynamic  para¬ 
meters  only.  In  the  linear  approach,  the  dynamic  parameters, 
such  as  the  moving  masses,  have  to  be  closely  matched  to 
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the  thermal  parameters  of  importance;  therefore,  optimization  of 
the  refrigerator  is  a  thermodynamic/dynamic  process,  more 
difficult,  of  course,  than  the  purely  thermodynamic  one.  Some 
aspects  of  the  thermodynamics  and  dynamics  of  the  proposed 
linear  approach  are  presented  in  Appendix  B. 

2.1.3  Passive  Vibration  Absorber 

The  use  of  a  passive  vibration  absorber  is  proposed  as  a  means 
for  eliminating  the  refrigerator  vibration  at  the  fundamental 
frequency.  A  vibration  damper  in  its  most  simple  form  consists 
of  a  mass  connected  to  a  spring  which  is  connected  to  a  machine 
for  which  attenuation  of  vibration  is  desired.  The  counter 
weight  of  a  linear-drive  Stirling  refrigerator  acts  as  a  bal¬ 
ancer  for  the  combined  sinusoidal  motion  of  the  piston  and 
displacer.  For  a  properly  balanced  system,  the  spring  force 
produced  by  the  reciprocating  counter-weight  will  be  just 
sufficient  to  nullify  the  forces  upon  the  machine  by  the 
moving  piston  and  displacer. 

Philips  has  successfully  applied  passive  vibration  absorbers  to 
cryogenic  refrigerators  and  demonstrated  a  reduction  of  cold 
finger  motion  by  nearly  two  orders  of  magnitude.  Specifically, 
the  MC-80  linearly  driven  refrigerator,  using  a  balancing  system 
similar  to  the  one  being  proposed,  has  a  cold  side  excursion  in 
the  range  of  3  to  5  pm.  It  is  probable  that  this  range  can  be 
lowered  with  more  sophisticated  balancing  techniques  than  those 
currently  used. 

2.1.4  Design  Parameters 

One  of  the  major  tasks  of  the  study  presented  in  this  report  was 
to  determine  the  thermal  and  physical  characteristics  of  a 
Stirling  refrigerator  design  that  can  meet  the  requirements 
associated  with  the  cooling  of  magnetometers  (and  other  super¬ 
conducting  elements).  Those  requirements  were  specified  in 
Section  C  of  the  ONR  RFP  No.  NOOOl  4-80-R-0007  as  follows; 
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.  Operation  at  a  temperature  of  10*K  (or  less) 
for  a  heat  load  of  50  mW. 

•  Electrical  input  power  not  to  exceed  250  W. 

•  Minimum  mechanical  motion  and  magnetic 
signatures . 

In  addition,  reference  was  made  to  8.5*K  as  the  desired  no-load 
temperature  and  to  4.2*K  as  a  longer  range  goal  for  the  oper¬ 
ating  temperature  of  small  cryogenic  refrigerators.  Also,  the 
variation  in  the  operating  temperature  level  was  to  be  main¬ 
tained  within  +  0.01 “K.  The  250  W  power  input  was  a  maximum, 
with  100  W  (or  less)  as  the  goal.  Finally,  the  refrigerator  was 
to  be  small,  light-weight,  and  reliable. 

To  meet  the  requirements  outlined  above.  Philips  Laboratories 
selected  a  triple-expansion,  linearly  driven,  balanced  design, 
as  discussed  in  the  previous  paragraphs.  Following  that  selec¬ 
tion,  an  optimization  study  was  conducted  to  arrive  at  re¬ 
frigerator  parameters  (i.e.,  working  pressure,  operating  fre¬ 
quency,  piston  and  displacer  strokes,  regenerator  dimensions, 
etc.)  which  would  yield  the  best  attainable  efficiency,  i.e., 
the  lowest  input  power  for  the  required  refrigeration. 

The  analytical  tools,  i.e.,  the  optimization  program  used  in 
the  study,  had  been  validated  through  years  of  experience  with 
operating  hardware.  However,  it  should  be  noted  that  designing 
a  cryogenic  refrigerator  for  very  small  cold  outputs,  such  as 
the  required  50  mW  at  lO’K,  is  an  extremely  sensitive  process: 
a  very  small  error,  or  unaccounted-for  loss,  can  wipe  out  the 
desired  net  refrigeration.  Therefore,  to  eliminate  that  danger, 
the  design  assumed  a  somewhat  higher  net  refrigeration  (0. 1  W  at 
10*K).  Thus,  the  results  of  the  analysis  are  conservative. 

The  major  refrigerator  parameters  that  resulted  from  the 
optimization  process  are: 
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Number  of  expansion  stages 

Mean  charge  pressure 

Working  gas 

Operating  speed 

Drive  system 

Type  of  linear  motor 

Piston  amplitude  and  dia. 

Displacer  amplitude 

Displacer  diameters 

Mass  of  piston  and  moving  coil 

Displacer  mass 

Displacer  spring  stiffness 

Linear  motor  efficiency 

Counter  weight 

Maximum  cold  production 

Maximum  power  input 

Cooler  weight 

Cooler  size 


5x10^  N/m^ 
helium 
SOO  cpm 

free  displacer  linear  drive 
moving  coil,  stationary  SmCo  magnets 
43.0  mm;  45.0  mm 
7.  0  mm 

14.3,  21.7,  31.6  mm 

3.0  kg 

0.3  kg 

8500  N/m 

70% 

passive  spring  mass  system 
0. 1  W  at  10®K 
200  W  electric 
74  kg 

630  mm  long,  350  mm  diam  (max) 


Two  items  are  to  be  noted.  The  parameters  listed  above  are 
subject  to  trade-offs;  for  instance,  efficiency  (or  power  input) 
can  be  traded  off  for  weight  (and  size).  Such  a  study  can  be 
conducted  in  the  next  phase  of  the  program  in  order  to  arrive  at 
an  "optimum  combination"  of  parameters,  rather  than  just  opti¬ 
mum  efficiency.  Also,  should  the  performance  evaluation  of 
the  hardware  indicate  that  the  design  was  conservative,  changes 
could  then  be  made  which  would  reduce  the  currently  predicted 
200  W  input  power  level  and  the  size  and  weight  of  the  cooler. 

2.1.5  Design  Layout 

The  data  that  resulted  from  the  analytical  studies  described 
in  the  previous  paragraph  was  translated  into  the  preliminary 
design  depicted  in  the  layout  drawing  of  Figure  3.  The  drawing 
contains  all  the  elements  of  the  proposed  refrigerator,  except 
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those  items  which  are  normally  designed  or  selected  during  the 
detailing  phase. 

The  major  features  of  the  cold  side  are  discussed  in  Par.  2.1.1. 
In  addition  to  those  features,  the  layout  drawing  addresses  the 
methods  for  measuring  the  thermal  performance  of  the  refrig¬ 
erator.  Specifically,  the  layout  shows  a  cavity  at  the  cold 
end,  which  is  to  serve  as  the  bulk  of  a  gas  thermometer,  the 
most  accurate  means  of  measuring  the  lowest  operating  tempera¬ 
ture.  The  capillary  shown  as  penetrating  the  cavity  wall  would 
connect  the  thermometer  bulb  to  an  external  pressure  indicator. 
For  redundancy,  provisions  were  made  to  incorporate  a  platinum 
resistance  thermometer  into  the  cold  region;  the  well  shown  next 
to  the  gas  thermometer  cavity  would  house  the  platinum  resistor. 

The  chamber  that  surrounds  the  cold  finger  is  a  vacuum  dewar  to 
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provide  means  of  maintaining  a  pressure  level  below  10  Torr 
in  that  region,  and  hence  to  thermally  insulate  the  cold  finger. 
It  should  be  noted  that  the  geometry  of  the  cold  finger/vacuum 
dewar  combination  can  be  readily  modified  to  accommodate  a 
magnetometer  or  another  superconducting  element. 

The  "neck"  that  couples  the  cold  finger  to  the  body  of  the 
refrigerator  is  the  ambient  heat  exchanger.  Tubulations  are 
provided  for  the  inlet  and  outlet  of  the  cooling  medium. 

The  piston/linear  motor  combination,  discussed  in  Par.  2.1.2,  is 
enclosed  in  a  hermetic  housing.  The  only  penetrations  through 
that  housing  are  the  helium  fill  valve  and  the  electrical 
feedthroughs  for  bringing  power  to  the  motor.  The  material  of 
the  housing  will  be  either  aluminum  or  titanium;  the  selection 
would  be  made  in  the  next  phase  of  the  program. 

The  proposed  refrigerator  is  630  mm  long,  and  the  diameter  at 
the  largest  cross-section  is  350  mm.  The  weight  breakdown  is 
as  follows: 


Housing  (assumed  aluminum) 

15.0 

kg 

Motor  (perm,  magnet  and  iron) 

35.0 

kg 

Piston  assembly 

3.0 

kg 

Counterweight 

15.0 

kg 

Displacer 

0.3 

kg 

Cold  finger 

1.2 

kg 

Heat  Exchanger 

4.5 

M 

Total 

74.0 

kg 

The  layout  drawing  of  Figure  3  would  provide  the  basis  for 
generating  the  detailed  design  during  the  next  phase  of  the 
program. 
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Critical  Areas 


Philips  Laboratories'  experience  with  the  triple-expansion 
Stirling  unit  which  was  designed,  built,  and  tested  in  1969-70 
indicates  that  the  performance  predicted  for  the  proposed  re¬ 
frigerator  can  be  attained.  However,  experience  with  the  same 
hardware,  as  well  as  with  other  cryogenic  refrigerators,  sug¬ 
gests  three  potential  problem  areas:  the  structural  integrity 
of  the  lead  matrix  in  the  lowest  temperature  regenerator; 
reliability,  especially  the  rapid  wear  of  conventional  piston 
seals;  and  the  need  to  maintain  temperature  stability  to 
+  O.OTk. 

These  problem  areas,  including  proposed  solutions,  are  discussed 
in  the  three  paragraphs  which  follow.  Another  problem  area, 
that  of  magnetic  siyr. tures,  is  discussed  in  Par.  2.3. 

2.2.1  Lead  Regenerator  Matrix 

The  proposed  design  is  based  on  the  assumption  that  the  lowest 
temperature  regenerator  matrix  will  contain  lead  spheres. 
Philips  Laboratories  used  lead  spheres  in  its  triple-expansion 
Stirling  refrigerator,  as  well  as  in  the  triple-expansion 
Vuilleumier  refrigerator  developed  under  contract  to  the  U.S. 
Air  Force  (Flight  Dynamics  Laboratory,  Wright-Patterson  Field). 
The  latter  unit  was  developed  for  possible  spaceborne  use; 
therefore,  life  and  reliability  were  critical  parameters. 
During  endurance  testing  of  the  Vuilleumier  refrigerator, 
difficulties  were  encountered  with  the  lead  sphere  matrix.  The 
Hughes  Aircraft  Company  encountered  similar  problems  during 
tests  with  a  triple-expansion  Vuilleumier  refrigerator  which 
also  contained  lead  spheres  in  the  lowest  temperature  regen¬ 
erator  matrix. 

Specifically,  the  lead  spheres  disintegrated  after  being  loaded 
into  the  housing  and  subjected  to  tests,  and  the  resulting  dust 
clogged  other  critical  passages.  Philips  Laboratories  conducted 
an  investigation  to  determine  the  cause  of  the  problem,  and 


concluded  that  the  failure  was  initiated  by  excessive  oxidation 
of  the  lead  spheres  (Hughes  Aircraft  reached  a  different  con¬ 
clusion).  A  subsequent  change  in  processing  and  packing  of  the 

lead  spheres  resulted  in  a  matrix  endurance  period  of  over  4000  ^ 

i 

hours.  i 

! 

Shortly  after  the  problem  was  identified,  Philips  Laboratories 
proposed  a  solution  to  the  lead  matrix  problem,  namely  the  use 
of  lead  mesh  rather  than  lead  spheres.  The  solution  was  imple¬ 
mented,  tested,  and  found  effective.  However,  since  lead  mesh 
is  a  "special"  material,  and,  unlike  spheres,  not  available 
commercially,  spheres  would  be  less  costly. 

The  endurance  evaluation  of  lead  spheres  is  continuing.  If  the 
results  of  the  evaluation  are  positive,  the  lowest  temperature 
regenerator  of  the  proposed  design  would  contain  lead  spheres. 

Otherwise,  Philips  Laboratories  would  recommend  the  use  of  mesh. 

Our  tests  have  shown  that  the  thermal  performances  of  spheres 
and  mesh  are  equivalent.  Thus,  the  performance  predictions  for 
the  proposed  design  apply  in  either  case. 

2.2.2  Piston  Seal 

The  sealing  effectiveness  and  endurance  characteristics  of  the 
piston  seal  are  critical  to  the  efficiency  and  reliability  of 
the  proposed  (as  well  as  any  other)  refrigerator  design.  Most 
commercially  available  cryogenic  refrigerators  use  expansion- 
type  (i.e.,  spring-loaded)  rings  to  seal  the  piston.  Those 
rings  have  acceptable  sealing  ability,  but  are  subject  to  rapid 
wear.  Other  refrigerator  designs,  notably  those  being  developed 
for  possible  spaceborne  applications,  propose  the  use  of  clear¬ 
ance-type  seals  which  rely  on  a  very  small  gap  between  the 
surfaces  to  be  sealed  to  prevent  gas  leakage  under  dynamic 
conditions.  Clearance  seals  are,  of  course,  not  subject  to  wear 
since  there  is  no  contact  between  the  adjoining  surfaces,  but 
their  sealing  effectiveness  is  lower  than  that  of  expansion 
seals. 
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Unlike  most  other  refrigerators,  the  piston  of  the  proposed 
design  is  not  subjected  to  lateral  thrust;  the  lack  of  such 
thrust  is  conducive  to  very  low  wear  rates,  even  when  con¬ 
ventional  seals  are  used.  For  the  application  at  hand,  however. 
Philips  Laboratories  proposes  to  use  a  novel  seal  design  having 
very  favorable  endurance  and  sealing  characteristics.  The 
performance  of  the  proposed  seal  was  verified  by  testing. 

The  rate  of  wear  of  a  piston  ring  depends  on  the  product  of  its 
surface  pressure  exerted  against  the  mating  cylinder  wall,  and 
on  its  linear  speed.  For  best  performance,  it  is  desirable  to 
have  the  sealing  ring  gas-pressure-balanced,  so  that  it  is  not 
pressed  against  the  wall  by  the  pressure  wave  variation.  Also, 
the  ring  should  be  free-floating  radially,  so  that  the  in¬ 
accuracies  in  the  sideways  guidance  of  the  piston  by  the  linear 
motor  are  not  transmitted  to  the  ring  and  thus  cause  it  to 
wear. 

The  proposed  piston  seal  design,  described  with  the  aid  of 
Figure  4,  incorporates  both  of  the  above  features.  A  reinforced 
Teflon  sleeve  is  fastened  to  the  free-floating  piston  body  and 


Figure  4.  Piston  seal  design. 


is  machined  to  a  sliding  fit  with  the  cylinder.  At  each  end  of 
this  body  are  hardened  steel  balls  backed  up  by  hardened  thrust 
washers.  The  balls  have  no  lubrication  and  carry  the  gas 
pressure  forces;  they  allow  the  thrust  washers  to  move  radially 
with  the  piston  without  transferring  a  side  force  to  the  ring 
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body.  A  phenolic  separator  plate  is  used  as  the  ball  spacer.  A 
seal  is  made  between  the  ring  body  I.D.  and  the  piston,  using  a 
silicon  rubber  0-ring  that  is  initially  compressed  more  than  the 
expected  side  motion  of  the  piston.  The  relatively  low  pressure 
wave  of  the  proposed  refrigerator  design  is  easily  handled  by 
this  soft  0-ring  material,  yet  minimal  side  force  is  transmitted 
to  the  floating  ring  body.  Since  the  ring  body  is  a  rigid 
member,  it  is  inherently  pressure  balanced. 

Because  of  this  pressure  balance  and  the  free-floating  features, 
ring  wear  is  minimal.  Calculations  show  that  this  type  of  ring 
will  wear  less  than  0.025  mm  radially  in  ten  years.  This 
minimal  wear  is  due  to  the  bearing  friction  induced  by  gas- 
pressure  thrust  and  to  the  0-ring  radial  load  imposed  by  the 
piston  motion. 

The  proposed  piston-ring  design  was  tested  in  a  double-expansion 
refrigerator  for  1000  hours.  At  the  end  of  this  period  there 
was  no  measurable  wear  or  visual  evidence  of  it  (i.e.,  wear 
particles).  Also,  there  was  no  phase  shift  in  the  refrigerator, 
an  indication  that  there  was  no  leakage  past  the  piston  seal. 


2.2.3  Output  Temperature  Fluctuation 


( 1 )  Introduction 

As  a  result  of  the  nearly  harmonic  time  variation  of  both  the 
pressure  and  volume  of  the  expansion  space,  the  cold  pro¬ 
duction  in  a  Stirling  refrigerator  fluctuates  about  a  mean  value 
with  time.  This  in  turn  gives  rise  to  a  fluctuating  temperature 
at  the  cold  end.  An  analysis  was  conducted  to  determine  what 
parameters  influence  this  fluctuation,  what  the  magnitude  of 
the  temperature  fluctuation  in  the  proposed  design  are,  and  how 
the  specified  +  O.Ol'K  tolerance  could  be  attained. 

( 2 )  Rate  of  Heat  Transfer  at  Wall  of  Expansion  Space 

The  helium  gas  in  the  expansion  space  is  subjected  to  cooling  or 
heating  during  the  thermodynamic  cycle.  The  rate  of  cooling  or 
heating  is  equal  to  the  expansion  or  compression  work  done  on 
the  gas  by  the  motion  of  the  displacer  and  piston.  The  wall  of 
the  expansion  space  is  in  turn  cooled  or  heated  by  the  helium 
gas.  Due  to  the  small  heat  capacity  and  high  turbulent  mixing 
of  the  helium  gas  in  the  expansion  space,  it  is  a  conservative 
approximation  (resulting  in  a  higher  temperature  fluctuation)  to 
equate  directly  the  rate  of  cooling  or  heating  of  the  wall  to 
the  work  done  by  the  displacer  and  the  piston.  This  results  in 
the  following  heat  transfer  rate  Q  at  the  wall: 


with  p  =  p^  (1  -  2e  cos  ( wt  -  0) ) 

and  ^o  ^  ^  ^ 


toV 

0  =  -  Pj^[l  -  2ecos(wt  -  e)][ - sinwt] 


This  can  be  written  as 
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If  the  first  and  second  term  of  the  left-hand  side  of  Equation 
(1)  are  denoted  by  and  AQ,  respectively.  Equation  (1)  can 
be  written  as 


Q  =  +  AQ 

where, Q  =  time-averaged  heat  transfer  rate  and  AQ  =  fluctuation 
in  the  heat  transfer  rate.  Since  e  <<  1,  aQ  can  be  approximat¬ 
ed  by ; 


AQ 


2 


sinwt 


(2 


The  effects  of  the  other  terms  in  AQ  can  be  included  by 
superposition  of  the  results  of  the  following  analysis. 

(3 )  Mathematical  Modeling  of  Thermal  Behavior  of  Wall 
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For  steady  state  operation,  the  cooling  load  on  the  outside 
surface  of  the  wall  is  and  the  heat  transfer  rate  on  the 
inside  surface  of  the  wall  is  +  aQ.  The  boundary- value 
problem  for  the  temperature  is: 


pc 


p  3t 


0  >  X  >  d 


m  1 

^  ~  A  ^^m  at  X  =  0 

,  3T  _  ^m 

^  9x  X  at  X  =  d. 


sinwt 


at  X  =  0 


{  3b  ) 


For  large  values  of  d ,  the  solution  to  this  problem  behaves 
like  that  of  a  semi- inf inite  body,  and  the  steady  periodic 
solution  is: 


U(x,t) 
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The  depth  at  which  the  temperature  fluctuation  is  damped  to 
within  5%  of  the  surface  amplitude  is,  from  Eq.  (4), 

’‘5%  =  -'•24  V! 


Therefore,  if  d  >  4.24 /a/u.  Equation  (4)  is  a  close  approximate 
solution  to  Equations  3a-3c. 

The  amplitude  of  fluctuation  of  temperature  T  at  x  =  d  can  be 
obtained  from  Equation  (4)  as: 


T  =  U 


p  V  w  I — 
o  fa 

2Ak 
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From  which. 
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Equation  (6)  gives  the  minimum  thickness  of  the  wall  through 
which  the  temperature  fluctuation  is  damped  to  T. 


0 
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It  can  be  observed  from  Equation  (5),  for  a  given  amplitude  of 
AQ,  i.e.,  Pn,uVQ/2,  T  decreases  with  higher  values  of  w  and  k 
and  lower  values  of  a.  An  effective  wall  material  for  the 
purpose  of  damping  the  temperature  fluctuation  should  have  high 
k  and  low  «. 


( 4 )  Design  Considerations 

To  produce  cold  with  a  temperature  fluctuation  not  to  exceed 
O.Ol'K,  therefore,  enough  thermal  inertia  has  to  be  built  in  the 
tip  of  the  cold  finger.  According  to  Equation  (5),  a  tip  made 
of  material  with  high  thermal  conductivity  k  and  low  diffusivity 
o  provides  better  temperature  stability.  At  temperatures 
around  10 “K,  lead  possesses  the  above  desirable  combination  of 
thermal  properties. 

The  wall  thickness  of  the  tip  of  the  cold  finger  required  to 
damp  temperature  fluctuations  down  to  specified  values  are  plot¬ 
ted  in  Figure  5.  The  total  cold  production  Q.pQ.j.;^L  1*94  W 
is  considered.  Since  the  net  cold  production  (total  cold 
production  minus  losses  through  the  cylinder  walls  and  dis¬ 
placer)  is  only  a  fraction  of  the  total  cold  production,  curves 
for  lower  values  of  cold  productions  are  also  plotted  in  Figure 
5.  The  corresponding  temperature  drops  across  the  wall  thick¬ 
ness  are  shown  in  Figure  6.  Figures  7  and  8  are  the  same  as 
Figures  5  and  6  with  a  larger  heat  transfer  area  A.  It  can  be 
seen  that  an  increase  in  A  not  only  significantly  reduces  the 
thickness  required  to  achieve  a  certain  temperature  stability, 
but  also  reduces  the  temperature  drop  across  the  wall  thickness. 

With  a  designed  cold  production  of  100  mW  (approximately  0^/20^ 
at  10*K,  a  temperature  stability  of  0.01'K  can  be  obtained  by 
using  a  lead  tip  of  1  cm  thickness  and  a  heat  transfer  area 
with  a  diameter  of  2.86  cm  (twice  the  diameter  of  the  expansion 
space  of  the  3rd  stage).  The  temperature  drop  across  this 
thickness  is  less  than  0.1*K. 
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Figure  5.  Wall  thickness  vs.  amplitude  of  temperature 
fluctuation  at  various  heat  transfer  rates  (heat  trans' 
fer  area  =  1.61  cm  ). 
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Figure  6.  AT  across  wall  vs.  amplitude  of  temperature 
fluctuation  at  varjous  heat  transfer  rates  (heat  trans 
fer  area  =  1.61  cm  ). 
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Fig\jre  7.  Wall  thickness  vs.  amplitude  of  temperature 
fluctuation  at  various  heat  transfer  rates  (heat  trans¬ 
fer  area  =  6.42  cm  ). 


34 


■^T  ACIioSiS  iv/* 


•oj  02  -03  ■<><?.  -05  ct 

/MPUtvDE  T^t^f^EKArueE  FtL'tTufl7/(?v,  T 


Figure  8.  AT  across  wall  vs.  amplitude  of  temperature 
fluctuation  at  various  heat  transfer  rates  (heat  trans¬ 
fer  area  =  6.42  cm  ). 
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In  order  to  increase  the  heat  transfer  area  A,  the  tip  of  a 
diameter  larger  than  that  of  the  expansion  space  can  be  used 
(Fig.  9).  To  ensure  that  the  heat-flow  process  makes  fuller  use 
of  the  geometric  cross-sectional  area  of  the  tip,  a  copper  cap 
(higher  conductivity  than  lead)  can  be  added  to  spread  out  the 
path  of  heat  flow  (Fig.  9b). 


CA>  WllMOuT  COPPER.  CAP _ (l>  ?  WtTH  foppgg  CAP 

Figure  9.  Schematic  of  heat-flow  path  in  tip 
of  cold  finger. 


( 5 )  Conclusions 

The  specified  +  0.01*K  output  temperature  tolerance  can  be 
attained  by  the  interposition  of  a  copper-sheathed  lead  mass 
between  the  expansion  space  of  the  refrigerator  and  the  item  to 
be  cooled.  Information  was  developed  under  this  contract  on  the 
relationship  of  the  attainable  temperature  fluctuation  levels  to 
other  important  parameters,  such  as  temperature  drop,  cold 
production,  and  mass.  The  final  selection  of  the  optimum 
geometry,  materials,  mass,  etc.,  of  the  method  for  attenuating 
temperature  fluctuation  can  best  be  done  during  the  next  phase 
of  the  program,  when  information  on  the  device  to  be  cooled  will 
be  available. 


2.3 


Magnetic  Effects  of  Mechanical  Motion 


The  coupling  of  magnetometers  (and  other  superconductive 
elements)  to  cryogenic  refrigerators  poses  a  difficult  problem. 
Magnetometers  are  to  be  used  in  field  applications  for  measuring 
very  low  frequency  fields  of  1  pT  or  less.  Cryogenic  re¬ 
frigerators  generate  the  required  low  temperatures  by  moving  two 
masses,  and  that  motion  induces  both  vibrations  at  the  cold  end 
and  strains  in  the  adjoining  walls.  The  vibrations  and  strains 
thus  induced  can  give  rise  to  magnetic  fields  higher  than  those 
that  need  to  be  detected.  One  of  the  most  difficult  challenges 
in  the  development  of  cryogenic  refrigerators  suitable  for  use 
in  conjunction  with  superconductors  will  be  to  arrive  at  designs 
which  will  minimize,  if  not  completely  eliminate,  the  magnetic 
fields  generated  at  the  superconductor/refrigerator  interface. 

One  of  the  factors  which  make  the  problem  difficult  is  the  lack 
of  experience  with  this  type  of  interface.  The  only  significant 
work,  to  our  knowledge,  was  done  at  the  National  Bureau  of 
Standards  by  J.E.  Zimmerman  and  his  colleagues.  However,  it  is 
obviously  desirable  to  have  an  inertially  balanced  refrigerator 
design  with  a  rigid  (minimal  strain)  cold  side  structure.  The 
balancing  approach  is  discussed  in  Par.  2.1.3.  The  question  of 
rigidity  needs  further  elaboration. 

Since  electrical  conductors  moving  in  the  earth's  magnetic  field 
generate  eddy  currents  and  unwanted  magnetic  fields,  all 
cold  side  structural  materials  should  be  non-conductors,  i.e., 
non-metals.  Indeed,  reinforced  plastics  have  been  used  in  the 
past  for  cold  finger  and  displacer  structures.  However,  the  use 
of  plastics  in  the  working  volume  of  a  cryogenic  refrigerator  is 
undesirable  since  outgassing  products  contaminate  the  working 
fluid  and  contribute  to  performance  deterioration.  The  proposed 
design  was  therefore  based  on  the  use  of  titanium  at  the  cold 
end  and  lead  and  phosphorous  bronze  for  the  regenerators. 
To  evaluate  the  field  that  would  be  generated  by  the  titanium 
moving  in  the  earth's  field,  a  preliminary  study,  very  general 
in  its  approach,  was  conducted. 
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As  indicated,  the  motion  of  the  electrically  conducting  com¬ 
ponents  of  the  proposed  refrigerator  in  the  earth's  magnetic 
field  will  produce  magnetic  interference.  The  conductors  react 
to  time-varying  magnetic  flux  and  to  accelerations  in  the 
presence  of  a  uniform  magnetic  field  by  generating  eddy  currents 
which  produce  reaction  magnetic  fields.  Worst-case  effects  are 
modeled  below,  and  estimates  of  magnetic  interference  are  made 
for  a  miniature  cryocooler  design. 

The  first  case  to  be  considered  is  that  of  a  pulsating  cy¬ 
lindrical  shell  (the  cold  finger,  for  instance)  in  a  uniform 
magnetic  field.  The  pulsating  shell  models  the  strain  of  the 
expansion  chamber  enclosure  induced  by  the  internal  pressure 
variations.  For  a  pressure  variation  of  the  form 

AP(t)  =  AP  sin  wt, 

a  thin  shell  will  respond  by  straining  an  amount: 

8r(t)  _  r  AP(t) 

r  ~  t  E 

where,  t  *  shell  thickness 
r  =  radius 

Sr  =  change  in  radius 
E  =  Young's  modulus 

Thus,  the  flux  area  normal  to  the  axis  varies  as: 

dA  _  2iTr^a)AP  cos  wt 
dt  “  tE  ' 

and  the  peak  magnetic  flux  change  is  thus: 


dt 


max 


2TTr  u)AP 
tE 


B. 
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where  B  is  the  magnitude  of  the  ambient  magnetic  field  in  the 
axial  direction.  The  changing  magnetic  flux  induces  a  cir¬ 
cumferentially  circulating  current  in  the  shell,  which  is 
inversely  proportioned  to  the  circuit  resistance  R.  Since  the 
reaction  magnetic  field  is  proportioned  to  the  circulating 
current , 

2 

u  B  r  0)  ,  T, 

=  o  ^  AP 

“reaction  ~  2p  E 


where  P  =  resistivity  of  shell 

Uq  =  magnetic  permeability  of  free  space. 


The  preceding  estimates  the  peak  field  at  one  end  of  an  axially 
semi- inf inite  solenoid  formed  by  a  thin  shell.  The  semi-infinite 
solenoid  approximation  is  conservative,  since  a  shorter  solenoid 
produces  a  weaker  field.  Also,  the  thin  shell  approximation  is 
somewhat  conservative,  since  the  shell  is  actually  constrained 
at  the  ends,  which  reduces  the  radial  strain  and  thus  reduces 
the  flux  area  variation. 


Parameters  for  estimating  this  effect  for  the  cold  finger  of  the 
proposed  refrigerator  are  as  follows: 


r  ~  1  cm 

AP  -  2  X  10^  N/m^ 

0)  ^  20  TTsec  ^ 

p  ~  5  X  10  ^  Q-m  (for  titanium  at  10 ’K) 

E  ~  1  X  10^^  N/m^  (for  titanium  at  10 *K) 


B  =  B  ~  50  |jT  for  worst  case  operational 

orientation  (horizontal  in  N-S 
direction ) 


Then, 


B 


reaction 


10  pT. 


A  field  of  this  magnitude  is  expected  to  present  a  problem  for 
measurements  in  the  pT  range.  However,  the  field  could  be 
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significantly  reduced  by:  shielding,  choice  of  operational 
orientation,  choice  of  material  with  different  p  and  E,  and 
variations  in  design  parameters  of  r,  t*)  and  AP. 

The  second  case  to  be  considered  is  that  of  a  changing  magnetic 
flux  due  to  motion  in  a  non-uniform  field.  In  unshielded  rooms, 
field  gradients  of  10  uT/M  are  not  uncommon.  This  effect  is 
estimated  for  the  displacer,  which  is  modelled  as  a  cylindrical 
shell  which  oscillates  in  the  axial  direction,  x,  with  amplitude 
a.  The  maximum  flux  variation  occurs  when  motion  is  along  the 
field  gradient.  The  cases  of  field  orientation  parallel  and 
normal  to  the  axis  of  symmetry  are  considered  separately  below. 

For  the  case  of  a  field  and  a  field  graaient  that  are  both 
axially  aligned,  the  flux  variation  for  the  previously  described 
cylindrical  shell  is: 


d(}) 

dt 


Trr 


2  ^ 
dt 


and 


dB 

at 


max 


jdB’ 'dx'  _  dB  ^  , 
|dX; ! dt!  dX 


Applying  the  semi- inf inite  solenoid  approximation 
yields  the  relation: 


B 


reaction 


dB 

4p  dX 


as  before 


For  the  same  parameters  as  before,  with  a  shell  thickness  t 
of  0.2  mm,  an  amplitude,  a,  of  6  mm,  and  a  field  gradient  of 
10  pT/M  in  the  axial  direction,  the  reaction  field  is: 

Bj.  ~  50  pT. 

The  effect  of  the  regenerator  matrix  within  the  oscillating 
shell  is  difficult  to  determine,  since  the  electrical  conduc¬ 
tivity  of  the  matrix  could  vary  drastically.  The  matrix  can  be 
modelled  as  a  set  of  concentric  shells  which  behave  as  the  shell 
analyzed  above.  The  reaction  field  is  then  proportional  to  the 
effective  matrix  conductivity.  The  conductivity  could  be  made 
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much  lower  than  that  of  the  outer  titanium  shell.  As  a  rough 
upper  bound,  the  regenerator  may  produce  interference  on  the 
same  order  as  that  of  the  outer  shell  described  above. 

For  the  case  of  axial  field  gradient  in  a  field  normal  to  the 
axis,  eddy  currents  will  circulate  about  the  field  direction  on 
the  surface  of  the  shell,  as  opposed  to  the  circumferential 
currents  of  the  first  two  cases.  The  cylinder  observes  the 
same  magnitude  of  time  varying  magnetic  field  as  the  prior  case, 
but  the  flux  area  is  now  the  projection  of  the  cylinder  surface 
area.  The  result  of  integrating  over  concentric  current  loops 
is  that  peak  reaction  fields  are  roughly  the  same  as  those  of 
the  prior  case  with  the  same  functional  dependence  on  r,  a, 
t,  P,  and  dB/dX. 

The  estimated  reaction  field  is  worst  case  since  it  assumes  that 
no  attempts  have  been  made  to  reduce  ambient  field  gradients, 
and  that  a  strong  gradient  is  aligned  with  the  machine's  axis. 
Also,  the  reaction  could  be  reduced  somewhat  by  varying  design 
parameters . 

The  last  effect  to  be  estimated  is  that  of  the  reaction  field 
due  to  the  acceleration  of  charges  in  a  uniform  magnetic  field. 
An  estimate  can  be  made  by  considering  a  flat  plate  oscillating 
in  the  presence  of  a  magnetic  field,  B^.  The  plate  consider¬ 
ed  oscillates  normal  to  its  surface.  The  magnetic  field  is 
directed  parallel  to  the  surface.  The  plate  has  dimensions 
by  "I  "  and  it  oscillates  with  amplitude  "a"  and  angular 
frequency  "u)".  Conduction  electrons  in  the  plate  experience 
forces  normal  to  the  velocity  and  to  the  field.  This  causes  a 
nearly  linear  charge  distribution  on  the  plate,  which  changes  as 
the  velocity  changes.  Redistribution  of  charges  results  in 
currents  in  the  plate  which  are  proportional  to  the  plate's 
acceleration.  The  resulting  reaction  magnetic  field  is  pro¬ 
portional  to  the  currents.  The  reaction  field  B^  is: 

Br  ~ 
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where,  =  permittivity  of  free  space  and  aw  =  peak 
acceleration  of  the  oscillating  plate.  Filling  in  constants  and 
rough  design  parameters  yields  an  immeasureable  response  of 
about  10  Bg,  or  about  10  T. 

The  preceding  analysis  suggests  that  the  influence  of  induced 
eddy  currents  could  be  a  problem.  It  appears  that  the  estimat¬ 
ed  interference  is  reducible,  possibly  to  acceptable  levels. 
This  preliminary  estimate  also  suggests  that  magnetic  noise 
could  result  from  the  operation  of  cryogenic  refrigerators  due 
to  vibrations  and  rotations  transmitted  to  the  magnetic  detector 
in  the  presence  of  the  earth's  field.  Thus,  the  use  of  con¬ 
ductors  as  refrigerator  components  should  not  be  ruled  out, 
since  structural  rigidity  and  vibration  compensation  are  likely 
to  be  of  critical  concern. 

It  should  be  noted  that  if  further  investigation  of  the  titanium 
cold  end  indicate  significant  problems,  the  titanium  can  be 
replaced  with  a  reinforced  plastic;  the  major  parameters  of  the 
design  would  not  change  significantly.  The  conduction  losses 
would  increase  somewhat,  and  the  unit  would  probably  have  to  be 
purged  periodically  to  "flush  out"  the  contaminants. 

It  should  also  be  noted  that  the  regenerator  matrices  now 
reciprocating  with  the  displacer  could  be  made  stationary.  In 
such  an  approach,  the  matrices  would  be  placed  outside  the 
displacer  in  an  annular  geometry.  Such  a  design  is  associated 
with  more  dead  space  and  with  a  more  complex  cold  side  construc¬ 
tion  . 

The  two  alternative  approaches  just  outlined,  i.e.,  the  use  of 
plastics  rather  than  titanium,  and  of  stationary  rather  than 
moving  regenerators,  are  practical,  but  thermally  and  struc¬ 
turally  less  attractive.  A  final  choice  can  best  be  made  as 
part  of  a  rigorous  trade-off  study  proposed  for  the  next  phase 
of  the  program. 
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2.4 


Remarks 


Following  completion  of  the  proposed  design,  an  effort  was 
made  to  compare  its  major  features  to  those  of  other  cryogenic 
refrigerators  which  operate  at  approximately  the  same  tempera¬ 
ture  level.  However,  the  only  refrigerators  which  operate  in 
the  8*K  to  15*K  range  are  very  large  units,  with  significantly 
higher  cold  producing  capacities  than  the  50  mW  output  of  the 
proposed  design.  Therefore,  the  comments  which  follow  are 
qualitative  rather  than  quantitative,  reflecting  the  difficulty 
of  a  one-to-one  comparison. 

The  proposed  design  is  based  on  the  Stirling  cycle,  which  is 
thermodynamically  reversible;  hence,  it  has  the  highest  attain¬ 
able  efficiency  (ideally).  Indeed,  experience  has  shown 
that  for  small  refrigeration  requirements,  no  other  cycle  has 
been  mechanized  with  higher  efficiency.  Our  study  suggests  that 
the  specified  efficiency,  i.e.,  50  mW  of  refrigeration  at  lO'K 
with  250  W  (maximum)  power  input,  can  be  readily  met.  Further¬ 
more,  the  100  W  ultimate  goal  appears  to  be  realistic. 

The  proposed  design  is  a  "one-piece"  configuration,  rather  than 
"split".  The  latter  approach,  whether  based  on  the  Stirling  or 
the  Gi f f ord-McMahon  cycle,  has  to  contend  with  significant 
losses  through  the  required  connection  between  the  expansion  and 
compression  spaces.  The  supposed  advantage  of  a  low-vibration 
expander  of  a  split  cycle  is  matched  by  the  very  low  vibrations 
attainable  with  a  balanced  one-piece  construction. 

A  significant  advantage  of  the  proposed  design  (versus  split 
cycle  approaches)  is  the  lack  of  valves,  which  are  associated 
with  losses  and  present  reliability  problems. 

It  should  be  noted  that  the  only  possible  disadvantage  of  the 
proposed  system  is  its  comparatively  high  weight,  which  is 
due  to  the  nature  of  a  slow-speed  linear  motor  and  the  need  for 
a  balancing  countermass. 
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3. 


OUTLINE  OF  PHASE  II 


3 . 1  Introduction 

Philips  Laboratories  recommends  that  the  design  presented  in 
this  report  be  detailed  during  the  next  phase  of  the  program. 
Those  details  would  then  be  translated  into  working  hardware 
during  a  subsequent  phase.  The  proposed  efforts  would  consist 
of  five  major  tasks: 

•  Parametric  trade-off  study  and  review 
of  current  design. 

•  Preparation  of  manufacturing  drawings. 

•  Review  of  material(s)  availability. 

•  Preparation  of  a  manufacturing  schedule 
and  a  test  plan. 

•  Report  of  Phase  II  results. 

The  suggested  tasks  are  discussed  in  detail  next. 

3.2  Discussion  of  Proposed  Program 

( 1 )  Parametric  Study  and  Review 

Prior  to  preparation  of  manufacturing  drawings,  efforts  are 
recommended  to: 

•  Refine  the  current  design  so  that  it 
reflects  the  Contracting  Agency's  priorities 
regarding  refrigerator  characteristics. 

As  indicated  in  this  report,  the  Philips  computer  program,  which 
mathematically  formulates  and  interrelates  the  various  design 
parameters  of  the  Stirling  cycle,  would  permit  the  selective 
optimization  of  the  cycle’s  important  characteristics.  For 
instance,  efficiency,  normally  an  optimized  parameter,  can  be 
traded-off  for  weight;  life  can  be  traded-off  for  weight;  etc. 
Philips  Laboratories  recommends  that  before  the  layout  drawing 
of  Figure  3  is  translated  into  details,  a  trade-off  study  be 
made  to  arrive  at  an  optimum  combination  of  important  re¬ 
frigerator  characteristics. 
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Such  a  study  rquires  a  priority  listing  of  refrigerator  char¬ 
acteristics,  i.e.,  a  ranking  of  efficiency,  weight,  size,  life, 
reliability,  etc.  Based  on  such  a  listing,  curves  can  be 
prepared  showing  the  interrelationships  of  the  various  param¬ 
eters.  With  the  assistance  and  concurrence  of  the  contracting 
agency,  nominal  design  points  would  be  selected  from  the  para¬ 
metric  curves. 

With  the  new  values  as  a  starting  point,  the  current  design 
would  be  reviewed  and,  most  likely,  minimal  dimensional  changes 
would  be  made.  These  changes  would  be  incorporated  in  the 
layout  drawing,  which  would  then  serve  to  initiate  the  next 
task,  viz .  , 

•  Make  final  choices  where  design  options 
still  exist. 

As  indicated  in  the  report,  many  design  choices  are  best  made 
during  the  transition  from  the  layout  drawing  to  the  manu¬ 
facturing  drawing.  Areas  which  would  require  further  consider¬ 
ation  during  the  next  phase  include;  the  lead  matrix  in  the 
coldest  regenerator,  the  method  of  static  sealing  (whether  with 
elastomeric  or  metal  seal),  and  the  material  for  the  refriger¬ 
ator  envelope.  Then, 

•  Expand  current  layout  drawing  to 
include  all  hardware  details. 

With  the  major  design  options  finalized,  the  next  task  of  the 
proposed  program  would  include:  selection  of  an  assembly 
procedure,  which  requires  that  fastening  methods  be  established; 
selection  of  motor  manufacturing  processes;  and  a  stress 
analysis  of  the  critical  refrigerator  components. 


( 2 )  Preparation  of  Manufacturing  Drawings 

Completion  of  the  task  outlined  above  would  permit  the  initia¬ 
tion  of  the  design  detailing  process.  It  should  be  noted 
that  the  detailing  process  includes  close  engineering  moni¬ 
toring  to  insure  that  the  thermodynamic  characteristics  of 


the  refrigerator  are  not  altered  by  design  compromises.  The  end 
result  of  this  task  would  be  a  drawing  package  from  which  an 
engineering  prototype  could  be  fabricated. 


( 3 )  Review  of  Material(s)  Availability 

There  is  a  possibility  that  some  of  the  materials  to  be  included 
in  the  proposed  refrigerator  are  not  available  "off-the-shelf". 
This  condition  normally  leads  to  significant  delays  in  the 
manufacturing  process.  To  eliminate  or  ease  that  condition, 
it  is  recommended  that  parallel  to  the  detailing  process, 
efforts  be  conducted  to  identify  long-lead  materials,  and  to 
place  purchase  orders  if  warranted.  This  activity  requires 
coordination  with  the  contracting  agency. 

( 4 )  Preparation  of  Manufacturing  Schedule  and  Test  Plan 

This  task  is  recommended  as  a  means  of  facilitating  the  planning 
and  conduct  of  the  subsequent  phase  of  the  program.  The  pre¬ 
paration  of  a  test  plan  can  be  especially  effective  in  this 
respect,  for  it  identifies  provisions  which  have  to  be  made  in 
the  hardware  for  accommodating  instrumentation.  Concurrently,  of 
course,  it  identifies  the  required  instrumentation  and  test 
facilities. 

( 5 )  Report 

The  fifth  task  would  be  the  preparation  of  a  report,  which  would 
include  a  description  and  discussion  of  the  design  developed 
during  this  phase,  and  detailed  plans  for  the  third  phase  of 
the  program. 
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Expansion  Staging  for  Regenerative  Cycles 
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APPENDIX  A 


EXPANSION  STAGING  FOR  REGENERATIVE  CYCLES 


Regenerative  type  refrigeration  cycles,  e.g. ,  Vuilleumier  (VM) 
and  Stirling,  produce  cold  by  the  reversible  expansion  of  an 
ideal  gas.  The  number  of  expansion  stages  and  their  geometry 
have  a  significant  effect  on  the  efficiency  with  which  the  re¬ 
quired  cold  can  be  produced.  The  design  aspects  of  single  and 
multiple  expansion  refrigerators  are  discussed  below. 


1.  Single  Expansion 

The  characteristics  of  the  VM  and  Stirling  cycles  are  normally 
described  by  depicting  one  expansion  space  of  the  refrigerator 
(Ref.  1) .  The  configuration  of  a  typical  VM  or  Stirling  cycle 
"cold  cylinder",  i.e.,  expansion  space,  displacer,  and  regen¬ 
erator,  is  shown  in  Figure  1.  In  principle,  ideal  regeneration 


DISPLACER 
GAS  FLOW 


Figure  1:  Typical 
Single  Expansion 
Configuration 


is  possible  in  both  cycles,  since  the  same  quantity  of  gas  flows 
through  the  regenerator  in  both  directions,  with  the  same  tem¬ 
perature  difference,  i.e.,  Tq  -  Tg.  The  amount  of  heat  rejected 
and  absorbed  by  the  gas  for  both  directions  of  flow  is  the  same, 
since  the  specific  heat  for  nearly  perfect  gases  is  essentially 
independent  of  pressure. 
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In  actual  practice,  the  regeneration  is  not  ideal,  i.e.,  there 
is  a  small  departure  from  the  ideal  caused  by  imperfect  heat 
transfer.  In  a  typical  regenerator,  a  quantity  of  heat  Q^-  must 
be  absorbed  and  rejected  during  each  cycle.  Due  to  various  im¬ 
perfections,  this  amount  is  reduced  to  Tl^Qr'  where  11^  is  the 
efficiency  of  the  regenerator.  This  means  of  course  that  not 
all  of  the  available  heat  is  transferred  to  the  regenerator. 

The  rest,  that  is  (l-n^)Qy,  is  carried  along  with  the  gas 
through  the  regenerator,  so  that  the  temperature  of  the  gas 
flowing  into  the  expansion  space  is  somewhat  above  that  desired. 
This  remainder,  which  we  will  designate  as  LQr,  constitutes  the 
regeneration  loss.  Since  this  loss  must  be  made  up  from  the  re¬ 
frigeration  produced  Qg,  it  should  be  compared  with  Qe-  Calcu¬ 
lations  show  that: 


AQr 

Qe 


=  Cj.(l 


Tlr) 


Tc  -  Te 


(1) 


The  constant  Cr  depends  mainly  on  the  compression  ratio  in  the 
refrigerator;  for  small  units,  its  value  is  approximately  10. 

For  example,  if  T^  =  300°K,  Tg  =  75®K,  and  the  regeneration  loss 
1-Tlj.  =  1%,  we  find: 


hQr 

Qe 


30% 


This  means  that  a  1%  loss  in  regeneration  involves  a  30%  loss  in 
refrigeration.  It  is  therefore  clear  that  for  efficient  opera¬ 
tion,  the  regeneration  loss  in  a  VM  or  Stirling  unit  has  to  be 
minimized. 

In  practice,  refrigerators  having  only  a  single  expansion  chamber 
are  referred  to  as  "single  expansion"  units. 


2.  Possible  Improvements 

Equation  (1)  for  the  relative  regeneration  loss  suggests  two  pos¬ 
sible  methods  for  improving  the  refrigeration  process:  decreasing 
the  value  of  C^.,  and  increasing  the  value  of 


C-. 
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The  way  to  decrease  Cr  would  be  to  increase  the  compression  ratio 
since  depends  mainly  on  this  ratio.  At  present,  no  practical 
configurations  exist  to  achieve  this  in  either  the  VM  or  Stirling 
cycles. 

The  value  of  1]^,  about  99%  in  typical  regenerative  refrigerators, 
can  be  increased  somewhat.  Regenerators  with  efficiencies  of 
99.5%  and  higher  have  been  made.  However,  since  a  regenerator  with 
such  high  efficiencies  normally  has  very  high  flow  losses,  lower 
refrigerator  efficiency  will  result. 


3.  Double  Expansion 

To  avoid  the  disadvantages  cited  above,  Philips,  in  1963,  modi¬ 
fied  the  conventional  single  expansion  system  (Ref.  2) .  The 
modification  consisted  of  thermally  staging  a  number  of  expan¬ 
sion  processes  in  as  many  expansion  spaces.  The  method  is  some¬ 
what  analogous  to  the  Keesom  cascade  process,  except  that  no 
separate  thermodynamic  cycles  are  used,  since  the  complete 
cycle  is  performed  in  one  closed  system.  This  modified  cycle 
has  been  referred  to  as  the  "double  expansion"  cycle. 

In  this  method  the  expansion  is  performed  in  two  expansion 
spaces,  each  at  a  different  temperature.  The  spaces  are  con¬ 
nected  in  series,  with  regenerators  in  the  interconnecting 
passages.  The  displacer  used  in  the  single  expansion  machine 
has  been  adapted  to  the  double  expansion  process,  as  shown  in 
Figure  2.  The  additional,  intermediate  expansion  space  is 
obtained  by  "stepping"  the  diameter  of  the  displacer.  By  vir¬ 
tue  of  this  stepped  displacer,  volume  variations  of  the  two  ex¬ 
pansion  spaces  (E  and  M)  are  in  phase.  The  advantage  of  this 
type  of  construction  is  that  the  pressure  difference  across  the 
seal  between  the  expansion  spaces  is  small,  as  in  the  conven¬ 
tional  design. 

As  a  result  of  the  additional  expansion  space,  the  refrigera¬ 
tor  produces  cold  at  two  different  temperatures,  whereas  in 
the  conventional  process  cold  production  takes  place  at  only 
one  temperature.  The  resultant  advantages  are  discussed  in 
Paragraph  4. 


The  following  discussion  shows  how  this  system  reduces  the 
influence  of  the  regeneration  loss.  Assume  that  a  mass  of 


COLO  SURFACE 


r 


EXPANSION 
SPACE  "E" 

REGENERATOR  2 

INTERMEDIATE 
EXPANSION  SPACE 

REGENERATOR  I 

SEAL 


M 


II 


Figure  2;  Typical  double  expansion  configuration 


working  gas  nij^  is  expanded  in  the  intermediate  expansion  space 
at  a  temperature  T^jj,  and  that  a  mass  of  working  gas  me  is  expanded 
in  the  top  expansion  space  at  a  temperature  Tg  -  both  gases  ex¬ 
panding  from  a  pressure  p^  to  a  pressure  P2 .  The  ideal  cold 
productions  in  the  expansion  spaces  are: 


The  mass  flowing  through  regenerator  1  between  room  temperature 
and  the  intermediate  temperature  now  is  m^  +  m^,  and  the  mass 
flowing  through  regenerator  2  between  the  intermediate  space 
and  the  expansion  space  is  mg.  The  regeneration  losses  for 
these  regenerators  are: 


I 


*°rl  =  l-V  +  ">£>  ■  ^rl>  ‘’’c  -  T„) 

‘^°r2  =  ”e<=p<1  -  \2'  (•'m  -  ’'e> 


Hence,  the  relative  regeneration  losses  are; 


°M  "m  \  ’ 


-^=  Crd  - 


(2) 


It  should  be  noted  that  Cj-  again  has  the  same  value  if  P1/P2 
has  the  same  value,  i.e.,  about  10. 

Now  AQjr2/QE  small,  because  is  of  the  order  of  140°K,  and 
Tm”T£  is  therefore  less  than  T^-Tg,  the  quantity  that  occurred 
in  Equation  (1).  The  consequence  is  that  little  of  the  cold 
production  in  the  top  expansion  space  is  lost,  thereby  result¬ 
ing  in  a  relatively  large  increase  in  net  cold  production. 

When  the  refrigerator  of  interest  needs  cold  at  one  level  only, 
no  net  production  Qm  is  needed,  and  it  is  permissible  for 
AQrl/QM  order  of  unity.  Taking  the  efficiency 

of  regenerator  1  to  be  98%,  it  follows  from  Equation  (2)  that; 


^Qrl  _  ^  ^E 


I 


Hence,  (m^j  +  m£)/m2^j  can  have  a  value  of  about  5.  This  means  that 
in  a  refrigerator  using  a  double  expansion  configuration  only  a 
mass  of  about  one-fifth  of  the  total  mass  flowing  through  the 
first  regenerator  needs  to  be  expanded  in  the  intermediate  space; 
this  causes  only  a  relatively  small  increase  in  the  compression 
work. 
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4. 


Applications 

It  is  evident  from  the  foregoing  that  given  the  same  conditions, 
i.e.,  cold  production  at  a  given  operating  temperature,  the  double 
expansion  configuration  is  more  efficient.  There  is,  however, 
additional  complexity,  viz. ,  two  regenerators  instead  of  one. 

It  should  be  noted  that  an  intermediate  temperature  level  in  the 
double  expansion  configuration  offers  significant  advantages  in 
some  cases.  For  instance,  in  certain  applications  where  radia¬ 
tion  shielding  is  desirable,  the  shield  can  be  cooled  at  the  in¬ 
termediate  temperature  level.  The  required  refrigeration  can  of 
course  be  produced  more  efficiently  at  this  level. 

Both  the  single  and  double  expansion  configurations  have  been  in¬ 
corporated  in  many  operational  systems.  Figure  3  shows,  for  com¬ 
parison,  single  and  double-expansion  displacer-regenerators  of  a 
VM  refrigerator  producing  1  W  at  77“K.  The  construction  is  simi¬ 
lar  in  both  cases. 


Figure  3:  Single  and  double 
expansion  displacer-regenera¬ 
tors  of  a  VM  refrigerator 
producing  1  W  at  77 °K. 


The  multi-expansion  concept  has  been  extended  one  step  further 
to  a  triple-expansion  Stirling  cycle  configuration  (Ref.  3),  as 
well  as  to  a  triple-expansion  Vuilleumier  refrigerator.  In  both 
cases,  temperatures  as  low  as  7.2*K,  with  only  one  cycle,  have 
been  attained  (Ref.  4).  Figure  4  shows  the  displacer-regenerators 
of  the  Vuilleumier  refrigerator. 
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Figure  4:  Displacer- 
regenerators  of 
Vuilleumier  refrigerator 
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APPENDIX  B 


Dynamics  of  a  Free-Displacer,  Free-Piston 
Stirling  Refrigerator 
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APPENDIX  B 


Dynamics  of  a  Fr ee-Di spl acer ,  Free-Piston 
Stirling  Refrigerator 


1 •  INTRODUCTION 

The  Stirling  cycle  is  characterized  by  the  interphased  reciprocat 
ing  motion  of  its  two  major  elements,  the  piston  and  the  dis¬ 
placer.  Traditionally,  this  motion  was  generated  by  either 
conventional  crank-type  mechanisms  or  by  the  Philips  rhombic 
drive.  In  both  instances,  the  rotary  motion  supplied  by  an 
electric  motor  had  to  be  translated  into  the  required  rectilinear 
motion  of  the  piston  and  the  displacer. 

In  the  early  1970 's,  the  Philips  Research  Laboratory  of  Eindhoven 
The  Netherlands,  designed  a  Stirling  cryogenic  refrigeration 
system  in  which  the  piston  was  reciprocated  by  a  linear  electric 
motor;  the  piston  in  turn,  by  fluidic  force  coupling,  drove 
(linearly)  the  displacer.  The  novel  approach  eliminated  the 
need  to  translate  rotary  to  linear  motion. 

The  new  refrigerator,  currently  in  industrial  and  military  use, 
offers  constructional  simplicity,  and  hence  high  reliability. 

Its  major  characteristics  are  described  in  the  following  section 
as  an  aid  in  understanding  the  design  presented  in  this  repo.t. 


2 .  DESCRIPTION 

A  simplified  cross  section  of  the  refrigerator  is  shown  in 
Figure  1.  The  low  temperature  (or  expansion)  side  contains  the 
displacer  which  contains  the  thermal  regenerator.  The  ambient- 
temperature  (or  compression)  side  contains  the  piston  which  is 
directly  coupled  to  the  moving  coil  of  the  permanent-magnet- type 
linear  lotor. 
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Fig.  1,  -  Schematic  cross 

section  of  linear 
refrigerator . 


The  permanent  magnet  maintains  an  approximately  constant 
field  in  the  gap  which  surrounds  the  inner  iron  yoke  of 
the  motor.  The  coil,  activated  by  an  AC  current,  recipro¬ 
cates  in  the  gap,  imparting  the  required  linear  motion  to 
the  piston.  A  spring  attached  to  the  outer  extremity  of  the 
piston/coil  subassembly  defines  the  center  position  of  the 
piston , 

The  displacer  is  driven  by  the  pressure  differential  estab¬ 
lished  by  the  losses  of  the  helium  working  gas  as  it  flows 
back  and  forth  through  the  regenerator  matrix.  A  spring 
attached  to  the  ambient  temperature  extremity  of  the  dis¬ 
placer  insures  proper  phase  relationship  between  the  moving 
masses . 

It  should  be  noted  that  all  the  forces  acting  on  the  piston 
and  displacer  are  in  the  direction  of  their  motion.  As 
a  result,  no  side  loads  are  imposed  on  the  piston  and  dis¬ 
placer  seals,  which  accounts  for  the  favorable  endurance 
characteristics  of  the  linear  Stirling  refrigerator. 
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3 •  PRESSURE  VARIATIONS  IN  A  STIRLING  REFRIGERATOR 

The  design,  and  hence  the  understanding,  of  this  type  of  Stirl¬ 
ing  refrigerator  requires  that  the  pressure  variations  which 
occur  in  its  working  volume  be  well  formulated.  A  theoretical 
expression  for  these  pressures  variations  is  given  in  References 
2  and  3,  and  is  not  discussed  in  this  report. 

The  pressure  variations  which  occur  in  the  refrigerator's 
working  volumes  are  induced  by: 

•  Motion  of  the  piston.  Its  position  determines  the  total 
volume . 

•  Motion  of  the  displacer.  Its  position  determines  the  amount 
of  gas,  at  a  given  time,  in  the  compression  (i.e.,  cold) 
space  and  the  expansion  (i.e.,  warm)  space. 

•  Gas  flow  induced  by  motion  of  the  piston.  This  flow  gives 
rise  to  pressure  differences  in  the  working  volumes. 

•  Gas  flow  induced  by  motion  of  the  displacer.  The  effects  of 
this  motion  on  the  pressure  in  the  refrigerator  are  similar 
in  nature  to  those  caused  by  the  motion  of  the  piston. 
Specifically: 

The  motion  of  the  piston,  with  the  displacer  in  a  fixed  position, 
induces  a  pressure  p  as  given  below: 

P  =  Pm  *  VO  y  =  Pm  ♦  ^ 

o  ' 


where,  p  =  average  cycle  pressure 

Ql 

S  =  piston  surface  area 

P 

y  =  piston  position 

=  maximum  expansion  volume 


V.  =  volume  of  element  i 
1 

T^  =  temperature  compression  space 
Tj^  =  temperature  of  element  i 
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with  the  piston  in  a  fixed  position  and  the  displacer 
reciprocating,  pressure  p  is  given  by: 


I>  S  (Z*-! 
m  n _ 


—  =  p  C 


(2) 


where,  S,  =  displacer  cross-sectional  area,  T  =  T  /T  the 
d  c  e 

temperature  ratio  between  the  two  varying  spaces,  and  x  = 
displacer  position^ 


The  change  in  pressure  caused  by  the  motion  of  the  displacer 
is  of  course  induced  by  the  working  gas  shuttling  between 
a  cold  and  a  warm  volume. 


The  pressure  in  the  cycle,  with  no  flow  losses  taken  into 
consideration,  is  given  by: 


P 


+  C  y 

y 


+ 


c 

X 


X 


(3) 


where:  x 


displacer  position 


y  =  position  position 

C  =  constant  related  to  pressure  variation 

induced  by  position  of  displacer 


C 

y 


constant  related  to  pressure  variation 
induced  by  position  of  piston. 


The  flow  losses  of  course  have  to  be  considered,  since  they 
affect  the  pressure  in  the  various  spaces  of  the  refrigerator. 
The  pressure  differences  across  a  resistive  element  (such 
as  a  regenerator  matrix  for  instance)  which  is  induced  by 
the  gas  flow  (assumed  to  be  laminar)  is  given  by 


Ap 


f'"'h  ‘’h  ^ 


c.v 


(4) 


where  : 

1  =  IpriRth  of  rnaistancc 

tljj  =  ch.irnc  t  I'rl  s  t  ic  dinmcter 
f  s  tloiisity 

V  =  ciixa  voloclty 
^  a  V  i  .s  c  o  s  i  f  y 
c  and  arc  constants* 
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It  should  be  noted  that  the  velocity  of  the  working  gas  at 
a  point  in  the  refrigerator  is  the  vectorial  sum  of  two 
velocities  which  are  proportional  to  the  velocity  y  of  the 
piston  and  the  velocity  x  of  the  displacer.  Also,  the 
pressure  difference  across  a  resistive  element  can  be  trans¬ 
lated  to  pressure  changes  on  its  two  sides. 

The  foregoing  analysis  and  discussion  leads  to  expressions 

giving  the  pressure  variations  (P^)  in  the  compression 

volume  and  that  (p  )  in  the  expansion  volume. 

e 


Pc  =  ‘’.n  *  S'-  * 


(5) 


Pe  Pm  *  ""yy  ^ 


(6) 


and  to  the  difference 


(7) 


where : 


a  constant  related  to  the  flow  loss  induced 
by  the  piston  velocity 


''d 

c  and  e 


a  constant  related  to  the  flow  loss  induced 
by  the  displacer  velcoity 


subscripts  referring  to  the  compression  and 
expansion  sides,  respectively 


PP 


C  -  C 
pe  pc 


'dd 


C  -  C 
de  dc 


X  =  displacer  velocity 
y  =  piston  velocity 


Equations  <5),  (6)  and  (7)  are  valid  only  if  constants  C  ,  c  , 

y 


C.,  and  C  are  independent  of  «,  y,  x,  and  y,  which  is  not 

dd  PP 


quite  the  case.  Consequently,  these  three  equations  are  approxi¬ 
mations  only.  As  shown  in  Figures  2,  3  and  4,  the  approximations 
are  acceptable. 


It  should  be  noted  that  in  addition  to  the  factors  which 
induce  pressure  variations  discussed  thus  far,  there  are 
several  other  second-order  contributing  factors, 
viz.,  leakage  past  the  piston,  temperature  variations 

due  to  heat  transfer  to  walls  and  due  to  imperfect  mixing, 
changes  in  the  average  temperatures  in  the  heat  exchangers 
and  in  the  varying  volumes  due  to  changes  in  loading,  and 
several  others.  In  addition,  the  flow  is  not  laminar  through¬ 
out  the  refrigerator.  Therefore,  flow  resistances  are  not 
linearly  dependent  on  gas  velocities,  as  assumed  in  Equation  (4) 


The  variation  of  constants  C 


C,,  and  C  was  investigated 
dd  pp 


with  the  aid  of  the  Philips  Stirling  Computer  Program  applied 
to  the  MMC-ao  refrigerator,  a  commercially  available  unit. 
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Figure  2  shows  the  variation 

of  C  and  C  as  a  function 
y  PP  ^ 

of  piston  amplitude  y  (with 

the  displacer  stationary) . 

As  the  graph  indicates, 

varies  only  slightly,  most 

likely  due  to  changes  .in 

temperature  at  different 

points  in  the  refrigerator. 

The  variation  in  C  is 

PP 

more  pronounced,  due  in 
large  part  to  non-laminar 
f  low . 
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Figure  3  shows  the  variation 
of  C  and  C,,  as  a  function 

X  dd 

of  displacer  amplitude  x 
(with  the  piston  stationary) . 
As  in  the  previous  case, 
the  variation  in  C  is  com- 

X 

paratively  small ,  and  that  of 
Cj,  is  more  pronounced. 
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Figures  4a  and  4b  show  the  variations  of  the  four  constants 
as  a  function  of  the  phase  angle  between  the  movements 
of  the  piston  and  the  displacer.  As  the  graphs  indicate, 
the  variations  are  relatively  small. 


4 .  COLD  PRODUCTION  AND  INPUT  POWER 

The  expressions  for  the  pressure  variations  which  occur  in 
the  refrigerator  serve  in  determining  its  cold  production 
and  the  power  input  to  it. 

The  ideal  cold  production  (Qg^)  is  given  by  the  expression: 

Ooo  =  n^p  dV^  =  -  nirCy  ».  f  f^r  (8) 

whore 

dVj  =  S^dx  =  Volviinc  variation  of  expansion 
volume 

n  =  number  of  cycles  per  second  ,  i.e.,  speed 
X  s  S  cos  ot 

y  -  y  cos  (ot.  <f>  ) 

ot  »  w  t 

U;  ■  angular  frequency 

As  Equation  (8)  indicates,  only  C^ ,  reflecting  the  pressure 
variations  induced  by  the  piston,  enters  into  the  expression 
for  the  ideal  cold  production  of  the  refrigerator. 


The  ideal  input  power  (N^)  is  given  by  the  expression; 


Np  =  n  ^  pdVp  =  n  *  y  Sp  sin  <f 


dr* 


(9 


where;  dv  =  s  dy, 
P  P 


As  can  be  seen  in  Equation  (9),  only  ,  or  the  pressure 
variations  induced  by  the  movement  of  the  displacer,  affect 
the  ideal  input  power  to  the  refrigerator. 

From  the  above  equation,  two  modes  of  operation  can  be  observ 
ed  for  a  positive  sin  <|>^^ ,  viz.,  shaft  power  is  put  into  or 
obtained  from  the  system  (C^  >  0,  <  0,  respectively). 

More  specif Ically ,  these  two  modes  are  the  classical  Stirl¬ 
ing  refrigerator  (t  >  1 )  and  Stirling  engine  ( t  <  1 )  cycles. 


5 •  MOTION  OF  FREE  DISPLACER 

The  most  significant  advantage  of  the  concept  depicted  in  Fig¬ 
ure  1  is  the  free  motion  of  the  displacer/  viz>,  no  mechanical 
means  are  required  to  reciprocate  it  in  the  desired  phase 
relationship  to  the  piston.  However,  to  insure  proper  motion, 
a  number  of  parameters  have  to  be  judiciously  established. 

Some  aspects  of  displacer  parameter  interaction  are  presented 
in  this  section. 

The  equation  of  motion  for  the  displacer  is  given  by; 

X  +  (p  -  p  )  S,  +  C  =  0  (10 

a  ecu  va 

where,  M^x  =  acceleration  force  (M.  is  the  displacer  mass), 
d  d 

(p  -p  )S  =  driving  force  on  the  displacer,  and  C  ,x  = 
e  c  d  yd 

spring  force  (friction  is  neglected). 


Integrating  the  expressions  for  the  pressure  into  the  above 
equation  yields: 


C  ,x 
vd 


a  (-C 

'  pp'  d 


(C  <0), 
'  PP 


(11) 


which  is  the  equation  for  a  damped  oscillator  with  a  driving 

force,  in  which  the  driving  force  (-  C  )  S^y  is  the  result 

PP  d-' 

of  the  piston  motion,  and  the  damping  force  C,,S,x  results 

dd  d 

from  the  motion  of  the  displacer. 

If  one  sets: 

y  =  y  cos  (a  -  and 

X  =  X  cos  a, 

then : 

where , 


COS'*  ♦  = 

-  9’ov'> 


(14) 


cor  = 


V<| 


(12) 

(13) 
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where 


i 


[ 

i 


If  0)  =  0)^,  which  IS  the  condition  where  the  displacer  is  tuned 
to  the  running  frequency,  =  0,  and  consequently  no 

refrigeration  is  produced.  However,  if  >  o) ,  then  (j)^^  >  0, 
and  the  unit  will  operate  as  a  refrigerator. 


yxo 


5.  -  Phnsc  shirt  anrl  rolatiV('  amplitude 
cf  the  displacer  versus  (U^/ uJ 


In  Figure  5,  x  and  ij)^^  are 
given  as  a  function  of 

As  the  graph  indi¬ 
cates,  when  01  /w  increases, 
w 

increases,  but  x  de¬ 
creases  . 

If  equation  (16)  is  inte¬ 
grated  into  (8)  , 


3  -nTrc  c  .s.?" 
eo  y  pd  d'^ 


sin  cos 


(17) 


and  the  optimum,  or  the  maximum  amount  of  refrigeration  will 

occur  at  d),  =  45®,  or 

^dr 


-nTTc  s  ,c 

_ V  <1  • 


(18) 


Then,  according  to  (16) 


Introducing  a  friction  force 
becomes 

M.  (w,'  -  CJ~) 

tan  O)  .  s  ; — 1 , ^ 

^dr  '•l/TT/S 


Fx/ I  X 1 ,  Equation  (15) 


(21) 


and  Equation  (16)  becomes 


*  =  Cpj?  cos  <f 


hylTT 


(22) 


With  the  introduction  of  the  friction  force,  the  amplitude 

X  decreases  slightly,  causing  a  corresponding  decrease  in 

Q  ,  which  will  now  occur  at  a  somewhat  smaller  value  of  d)  .  . 
eo  dr 


6 .  MOTION  OF  THE  PISTON 

As  indicated  in  Section  1,  the  piston  is  driven  by  a  linear 
electric  motor.  Its  equation  of  motion  is  given  by: 


V  * 


C  y 

vp^ 


(23) 


where 

M  y  is  the  acceleration  force  (M 
is  mass  nr^movin/'  part), 

force?  on  the  pi&ton,  C  y  the  spring  force 
find  D  I  i  the  e  1  ec  t  roma};no  t i c  force.  B  Is 
tlio  mflgnctic  field  in  ttie  gap,  I  the  ?ongth 
of  wire  present  in  the  gap  and  i  the  current 
vlirough  the  coil  DJ 


It  should  be  noted  that  the  friction  force  was  neglected  in 
equation  (23). 


The  electric  current  through  the  motor  coils  is  given  by 


U  cos  {/!*(/>)  .  u  U„co»(/»*/)  _  D  1 

1  *  -2 -  - m  ^  _2 - ^  j24) 

i  Ri 

whe  re 

s  the  voltngc  of  the  source 
and  ^ ia  the  faso  between  vol- 
on  movemoii t ) ,  U  is  the  induct- 
ivo  voltage  due  to  the  movement  of  the  coil 
ami  thus  equal  to  y. 


U  cos(  fl  *  V)  i 
(%.«-  ^ 
tage  and  pisl 
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When  Equation  (24)  and  the  expression  for  the  working  pressure 
are  integrated  in  (23), 

(M  W-.S  C  -Cv  )?  cos  /3  +  sin/3  = 

(25) 

SpC^.«  cos  (/3*  cos 


It  should  be  noted  that  the  above  expression  does  not  include 
the  flow  losses.  Also,  the  expression  includes: 

y  s  ^  COS  /3 
X  =  S  cos  (  /3  + 

(H  1^,)^  D 

and  D  =  - 1*  —  and  B  =  — -  have  been 

C  H .  U  H . 

1  X 

introduced . 


The  spring  constant  for  the  piston  has  two  components:  ^vp' 
the  mechanical  spring  constant,  which  is  relatively  small; 
and  Cy  Sp,  the  constant  for  the  gas  spring.  Solving  equation 


(25)  with  the  aid  of  (21)  and  (22)  gives 


tan  y  = 


p.  U>*  S  C  C  .sin  40,  cos  to 
c  n  X  T>d  <ir 

see  ,cos'<r  .  -  (M  u>"-C  S  rC  ) 

p  X  pd  y  dr  '  p  y  p  vp' 


As  the  above  expression  indicates,  (j)  is  independent  of  the 
applied  voltage  ,  but  strongly  dependent  on  resonant 
condition.  Also, 


n  l-  sin 

II  <) 


■■  ^,'\‘p..  -'•‘■'y’.ir  (27) 

which  indicates  that  the  amplitude  of  the  piston  motion  is 
proportional  to  the  applied  voltage. 

The  power  delivered  to  the  piston  is  given  by: 


'■  p  "  J  '  ("u"o'‘’"  'f 


The  electric  power  input  is: 


,  .1?,.  •‘■'i'’ 


I 
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» 


I 


ft 


» 


I 


> 


>. 


Hence, the  efficiency  is  given  by: 


X 

-H.-. 


A  n 


cl  (  1  ♦  a)  (  A  +  cos 


wi  th 


-Jl. 


The  maximum  efficiency  is  reached  at  ()>  =  90°. 


(30) 


For  optimum  and  that  is  for  lowest  power  input  for 

highest  refrigeration. 


c  s  +c 


M  =- 
p  . 


1* 


C,./2 


+  C 


a>  ‘ 


vp 


c  c  ,c 

p  s 

"■  P  Pc  ^Pm 
Co^ 


(31) 


Thus,  if  C  is  neglected,  the  mass  of  the  piston-coil 
^  2 

assembly  is  inversely  proportional  to  to  .  The  speed  of  the 

refrigerator  is  therefore  limited  by  the  piston/coil  mass, 

which  is,  of  course,  determined  by  design  or  constructional 

considerations.  The  only  manner  in  which  to  can  be  increased 

is  by  increasing  p  S  . 

m  p 


7 .  EXPERIMENTAL  VERIFICATION 

The  analysis  outlined  thus  far  was  verified  with  the  aid  of  the 
MMC-80  refrigerator,  after  inclusion  of  some  of  the  effects 
previously  mentioned  in  Section  3. 

In  Figures  6a,  6b  and  6c,  the  displacer  amplitude  x,  the  phase 

shift  ,  and  the  cold  production  Q  have  been  plotted 

Q  r  0 

against  to^/os.  In  the  tests  used  to  derive  the  plotted  data, 

the  piston  amplitude  y,  the  cold-end  temperature  T  ,  and  the 

e 

average  pressure  p  were  Itept  constant. 


t: 


6 
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Fi/7*  -  I'onip.iri  ^on  h.' i  r  1  .-n  1  .i  r  *‘<1  niicl 

o  x;>.' ri  in.*n  I  .-1 1  rrsult-.  Di'.pian-r 
.inipi  1  tuili- ,  r.i-'i'  -iliii't  .mil  cold 

production  voisu?-  Oj  J  uj 


When  a  friction  loss  of  F  =  0.5  N  is  assumed,  a  realistic  value 
in  this  case,  the  agreement  between  theory  and  practice  for  x 
and  ‘1’^^  are  close.  The  measured  cold  production  is  somewhat 
lower  than  that  predicted,  which  is  often  the  case  in  refriger¬ 
ators  in  which  the  value  of  the  heat  losses  is  comparable  to 
that  of  the  gross  cold  production. 

In  Figures  7a,  7b  and  7c,  the  electrical  input  power  N  ,  the 

e  1 

displacer  amplitude  x,  and  the  cold  production  Q  are  plotted 

e 

against  a  varying  average  working  pressure  P  ,  with  the  piston 

in 

amplitude  y  kept  constant.  It  should  be  noted  that  since  the 
electric  input  power  is  given  by: 


'  net*-?)  Vs 

^  ^  sin  f 


the  above  expression  will  show  a  minimum  at  the  point  where  the 
piston  is  in  resonance.  This  is  confirmed  by  Figure  7a. 
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KOMENCLATURE 

p  pressure 
T  temperature 
T  temperature  ratio 
S  surface  area 
M  mass 
V  volume 

IJ  total  volune,  normal]. sod  to  T  and  re¬ 
duced  ^-o  V 

j .  .  n  .  . 

X  displacer  position 
y  piston  position 
an^^ul.ir  frorpioncy 

if  faso  shift  >)etwe<Mi  piston  and  applied 
VO  1 1  af'o 

^dr  betwefii  piston  and  «lispl  acer 

n  cycles  per  soronct 


t 

ot 

Q 

N 

"o 


R. 

1 

B 

c 

B 

u 

A 

F 


time 
OJ  t 
v  ~  If 


dr 


fr 


PP 

"dd 

“pd 


1/t 

Cold  production 

power 

officio ncy 

amplitude  of  the  applied  voltage 
inductive  voltage  duo  to  the  movonirnt 
of  the  coil 

resistance  of  the  coil 
magnet  field  in  the  gap 
length  of  wire  present  in  the  gap 
electric  current 

constant  =  (n  1  )  /ll. 

''  G  v'  '  i 

constant  =  (n  1  )  /ll. 

e  v'  '  1 

ratio  between  N  and  B  Uo^)  /2 

friction  force  acting  on  the  displacer 
constant  of  pressure  variation  due  to 
displacer  position 

constant  of  pressure  variation  due  to 
piston  position 

constant  of  flow  loss  due  to  piston  ve¬ 
locity 

constant  of  flow  loss  due  to  displacer 
velocity 

r  C  -  C 
pe  pc 


=  C 


-  C 


■'de  ''dc 

*  -C  /C  .  , 
pp'  dd 

spring  constant 


Snbscri nt  s 

m  mean  value 

c  compression  side 
o  expansion  side 

i  volume  i 

d  displacer 

p  piston 

O  reference  value 
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